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■ Abstract 

We study the dynamics of a chain of coupled particles subjected to 
a restoring force (Klein-Gordon lattice) in the cases of either periodic or 
' Dirichlet boundary conditions. Precisely, we prove that, when the ini- 

I tial data are of small amplitude and have long wavelength, the main part 

■ of the solution is interpolated by a solution of the nonlinear Schrodinger 

equation, which in turn has the property that its Fourier coefficients de- 
cay exponentially. The first order correction to the solution has Fourier 
coefficients that decay exponentially in the periodic case, but only as a 

I ' power in the Dirichlet case. In particular our result allows one to explain 

, the numerical computations of the paper [BMP07| . 

(N ; 

^ ! 1 Introduction. 

00 ' 

f-*) I The dynamics of chains of coupled particles has been the object of a huge 

. number of studies, but only recently some numerical works (see jBGOSilBMPO?! ) 

. 5^ I have shown that the boundary conditions have some relevance on FPU type 

^ . investigations. The goal of the present paper is to study analytically the effects 

?— ( ' of the boundary conditions (BC) on the dynamics of a simple 1-dimensional 

model, namely the so called Klein Gordon lattice (coupled particles subjected 
to an on site restoring force). Precisely, we concentrate on the cases of periodic 
and of Dirichlet boundary conditions, and use the methods of normal form to 
study the dynamics. This leads to a quite clear understanding of the role of the 
boundary conditions and to an explanation of the numerical results of |BMP07| . 
On the contrary, our theory does not allow one to explain the results of |BG08| . 

More precisely, we study the dynamics of a large lattice corresponding to 
small amplitude initial data with long wavelength; we show that if the size A'^ 
is large enough and the amplitude e of the initial excitation is of order jj, 
then the solution z has the form 

z = fiz''{t)+ n^zi{t) , (1.1) 
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up to times \t\ < 0{p~^). In (|l.ip z° is interpolated by a solution of the non- 
linear Schrodinger equation (NLS) and has a behaviour which is independent of 
the BC. On the contrary Zi depends on the BC. Precisely, its Fourier coefficients 
decrease exponentially in the periodic case, but only as in the Dirichlct 

case. 

The theory we develop in order to give the representation provides a 
clear interpretation of the phenomenon. Indeed, it turns out that the normal 
form of the system is independent of the BC (and coincides with the NLS), 
but the coordinate transformation introducing the normal form has properties 
which are different in the periodic and in the Dirichlet case. In particular, in the 
Dirichlet case it maps sequences which decay fast into sequences which decay 
as This introduces the slow decay in the Dirichlet case. 

It should be pointed out that our result still depends on the size N of the 
lattice0 Nevertheless, we think that (within the range of validity of our re- 
sult) we clearly show the role of the boundary conditions and provide a good 
interpretation of the numerical results. 

The present situation has many similarities with the one occurring in the the- 
ory of the Navicr Stokes equation (see e.g. |Tem91| ). where it is well known that 
the spectrum of the solution depends on the boundary conditions. Moreover, we 
recall that a power law decay of localized object has been previously observed 
in nonlinear lattice dynamics in |DP03[ |Pey04[ IFla98[ IGF05| and that the con- 
nection between the nonlinear Schrodinger equation and the dynamics of long 
chains of particles was studied in many papers (see e.g. [Kal89| [KSM92| [Sch98| ) . 

The paper is organized as follows: in sect. [2]wc present the model, state our 
main result and discuss its relation with numerical computations. In sect. [3]we 
give the proof of the normal form construction. In sect. [4] we give the proof of 
the decomposition p.ip . Some technical details are deferred to the appendix. 
Each section is split into several subsections. 

Acknowledgments. This work was partially supported by MIUR under the 
project COFIN2005 "Sistemi dinamici nonlincari c applicazioni fisichc" . We 
thank Antonio Ponno, Luigi Galgani and Simone Paleari for many discussion 
and suggestions during the preparation of this paper. We also thank Sergej 
Flach and Michel Peyrard for some interesting discussion. 



^ However our normal form holds in the region of the phase points with small energy density 
e^, independently of N. The limitation e ~ comes from the fact that we are only able to 
study the dynamics of the NLS in this situation. 
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2 Main result 



In this chapter wc present the model, we recall some numerical simulations (see 
|BMP07j ) which clearly show the dependence of the metastable Fourier decay 
on the boundary conditions and we finally state the main theoretical result and 
use it to explain the numerics. 

2.1 The model. 

We consider a chain of particles described by the Hamiltonian function 

= Ef+E^fe)+E^fe-*-i) (2-1) 

where j runs from to in the case of Dirichlct boundary conditions (DBC), 
namely qo = qN+i ~ 0, while it runs from — {N+ 1) to N in the case of periodic 
boundary conditions (PBC), i.e. q-M-i — qN+i- The corresponding Hamilton 
equations are 

9, = -V\q,) - W'{q, - <Z,_i) + W'{q,+^ - q,) . (2.2) 
We recall that the standard Fermi Pasta Ulam model is obtained by taking 

2 3 4 

y = and W{x) = \ + ot\ + l3\- Here instead we will take 

V[x) = + ^ax^ + ^[3x'^ , W{x) = ^ax^ , a, /?, a > . (2.3) 
Explicitly our Hamiltonian has the form 

H = H0 + H1+H2, (2.4) 

j 

H^{p,q) := "Ey' ^2(p,g):=/3Ex ■ (2-6) 

J 3 

Remark 2.1. In the case where 

V{x) = V{-x) (2.7) 

the equations (j2.2[) with PBC are invariant under the involution qj 1— > —q-j, 
Pj l—^ —p-j. As a consequence the submanifold of the periodic sequences which 
are also skew-symmetric, is invariant under the dynamics. For this reason when 
(|2.7p is fulfilled the case of DBC is just a subcase of the case of PBC. This 
happens in the standard FPU model and also in the case of the Hamiltonian 
(|2.ip with the potential (|2.3p and a = 0. The case a 7^ is the simplest one 
where a difference between DBC and PBC is possible. 
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Consider the vectors 




fc= 1,...,7V, 

k = —^ — /V 

e,(j)= <; V^+i^-^V^+iy ' (2.8) 

k = 0, 



V2Ar+2' 



A: = -A^- 1, 



then the Fourier basis is formed by Cfc, fc = 1, . . . , and 6pD = in the case 
of DBC, and by Cfc, k = —N — 1, . . . ,N and 6pD = 1 in the case of PBC. Here 
we will treat in a unified way both the cases of DBC and PBC, thus we will not 
specify the set where the indexes j and k vary. Introducing the rescaled Fourier 
variables {pk , (jk ) defined by 

Pj = y2 V^PkSkij), qj = V -^ek{j), (2.9) 



fe 

where the frequencies are defined by 



2 / kn 



a.. = ^/l + 4asin^(^^), (2.10) 



the Hamiltonian Hq is changed to 

Ho = j:^k^^. (2.11) 



2.2 The phenomenon and its numerical evidence. 

Let us define the energy of a normal mode and its time average by 



1 f* 

- / Ek{s)ds . 
* Jo 



In the case of PBC the oscillators of index k and —k are in resonance, so the 
relevant quantity to be observed is the average {E^) = ^{{Ek) + {E^k)). 

Take an initial datum with all the energy concentrated on the first Fourier 
mode with energy density £ = Hq/N = 0.001. Integrating the system numeri- 
cally one can see that after a short transient the time averages of the harmonic 
energies relax to well defined steady values, which persist for very long times. 
In figures [1] we plot in a semi-log scale the time-average energies (Ek) (T) (or 
{Ek){t)) at time T = 10^ (subsequent to the relaxation time) as a function 
of the index k. The parameters in the two panels are a = 0.25, (3 = and 
a = 0, P = 0.25 respectively. In both the cases a = 0.5. In each distribution 
the dots refer to the DBC case while the crosses pertain the PBC onc[l. 



^Actually we plot only modes with odd index, since, as shown in | BMP07| the dy namics 
involves only modes with odd index t 
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Figure 1: Averaged harmonic energies distribution. DBC (dots) and PBC 
(crosses) with N = 511, a = .5, 5 = 0.001, T = 10^ Panel (a): a model 
with a = 0.25. Panel (b): /3 model with /? = 0.25. 



While in panel (b) one clearly observes a perfect overlapping of the expo- 
nential part of the decays, in panel (a) a sharp difference arises. Indeed, while 
the PBC solution is once more characterized by an exponential distribution, in 
the case of DBC one sees a richer behavior: at an energy approximately equal 
to 10~® there is crossover and a new regime appears. Nevertheless, a striking 
similarity among the exponential part of the two dynamics is evident. 

To describe more carefully the situation in the case of DBC we plot in figures 
[2] four different distributions of the quantities (Ek)-, in a semi- log and in a log- 
log scale respectively. They correspond to different values of the energy density 
(see the caption). In the first panel we plot the first part of the distribution: we 
notice that by decreasing the energy density, the slope of the exponential decay 
of the low frequencies increases. In the second figure, instead, we focus our 
attention on the second part of the distribution: we see that the corresponding 
curves are parallel. So a change of energy only induces a translation. Except for 
the last part, that we will interpret as due to discreteness effects, the curves are 
very well interpolated by a straight line giving a power decay with an exponent 
close to —6. A similar behavior is also obtained if one excites a few modes of 
large wave length. 

2.3 Explanation of the phenomenon. 

In order to state our main result we need a topology in the phase space. 
Definition 2.1. Let us define the spaces ^ of the sequences p = {pk} s.t. 

Mia ■■= A^Et^]''^'""' \P^\^ < ^ ' W ■■= l^l> ' (2.12) 

fe 
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Figure 2: DBC with parameters N = 511, a = 0.5, a 0.1, T 10^ 
Panel (a): distribution of (Ek) in semi-log scale. Energy densities: 
£ = 0.05,0.025,0.01,0.005,0.001. Panel (b): distribution of (Ek) in log-log 
scale. Energy densities: £ = 0.01,0.001,0.0001,0.00001. 

and the phase spaces Vsm ^1,a ^ ^l.a 9 (Pi 9)- 

The main part of the solution will be described by NLS, so we consider a a 
smooth solution lp{x^t) of the nonlinear Schrodinger equation 

^idrV ^ ~dx.,^ + -1^\^\^ , •= ^ (-^ - ' a;eT:=M/27rZ (2.13) 

For fixed r we will measure the size of a function by the norm 

k 

where ipk are the Fourier coefficients of (p, which are defined by 

p{x) ■.= Y,^kem ■ (2.15) 

fc 

here e^(x) is the continuous Fourier basis, 

cos kx k > 

^ = 0- (2-16) 

^ sin(-fca;) fc < 

Remark 2.2. The dynamics of is wcU known [FT87[ [GK03] . Precisely, if 

7 < 7* then Vcr > there exists < cr' < a such that from |j(^(x,0)||g ^ = 1 it 
follows ||iy9'^(a;, i)||g ^, < Cs ~ 1 for all times. 

Corresponding to ip{x,T) we define an approximate solution of the original 
model by 

z]{t) EE {p'^it),q^it)) := (Re (e^V(w, afih)) , Im (e^V(w, afih))) (2.17) 
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Our main result concerns the comparison between z°'{t) and the solution z{t) 
of the original system with initial datum 



z(0) fiz'^iO) 

Theorem 2.1. Assume \\tp^{x)\\^^ < 1, then\/T > there exists /i* > with 
the following properties: if < jj.^ then there exists ziit) defined for \t\ < T/i^^ 
such that 

zit) ^ fiz^it) + fi^ziit) (2.18) 
where is the approximate solution just defined and, 

\\Mt)\\s..' <C (2.19) 

with 



i < s , cr' > if PBC 



s < I , cr' = if DBC 



(2.20) 



The above result gives an upper estimate of the error zi {t) . We want now to 
compute it, at least for short times. To this end we assume, for simplicity, that 
i{P is purely imaginary, i.e. that the initial datum has zero velocity. Consider 
now the complex function ■(/;io with Fourier coefficients given by 

(V'io)fc [4e''^'=* - 3e2« - e-^" _ 6ie'"^* + 6i] (2.21) 

6v 2 

where <^{x) -~\Lp^{x)^, and let zio = (pio7'7io) be the sequence with 

Pw,j + 1910, 



V2 



"010 (w) 



Wc have the following 

Theorem 2.2. For any < & < 1, one has 

\\zl{t)-z,o{t)\\,^,,<C^i'-\ \t\<-^ (2.22) 

M 

The above theorems provide the interpretation for the numerical results of 
the previous paragraph. Referring to FIG. 1, we identify the exponential part 
of the distribution as due to the main part of the solution, namely z°, which is 
the same for both the boundary conditions. In the case of PBC, the inequality 
(|2.19p implies that also the error zi is exponentially decreasing; thus the whole 
solution is in particular exponentially localized in Fourier space. 

In the case of DBC the situation is different. Indeed the correction, namely 
zi, is ensured to have coefficients such that the series with general term |fcp*|zi_fe(t)|^ 
is convergent; this is very close to say that 



U2S+1 



7 



which, taking s very close to 5/2 essentially gives a power law decay like \k\ ^. 
Then theorem 12.21 shows that this is actually optimal, as seen by taking 



ip°{x)=ismx (2.23) 
(as in the numerical computations). Indeed in such a case one has 

c|, = (sin2x)2~l (2.24) 

so that (|2.22p shows that after a time of order 1 the energy of the k-th mode is 
of order /i^/Zc^ as shown by the numerics. 



3 The normal form construction. 

In this part wc introduce and use the methods of the normal form theory for the 
proof of our main result. Accordingly one looks for a canonical transformation 
putting the system in a simpler form. 

3.1 Preliminaries and main claim. 

We first need to introduce some notations: 

• z will denote a phase point. In particular a phase point can be represented 
using the coordinates {pj, qj) of the lattice's particles or the Fourier coor- 
dinates {pk,qk)- 

• In the phase space we will also use coordinates ij} defined by 
and, in real space 

■=^^kek{j) 

k 

• Given a Hamiltonian function H, wc will denote by Xh the correspond- 
ing Hamiltonian vector field. Thus if one uses for example the variables 
{pj,qj), one has 

Correspondingly we will write the Hamilton equations of a Hamiltonian 
function H by 

z = Xh{z) . 
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• The Lie transform generated by a Hamiltonian function x is the time 
one flow of the corresponding Hamilton equations, namely 

• The Poisson bracket {/, g} of two function /, g is defined by 

The normalizing transformation will be constructed by composing two Lie 
transforms 71 = and T2 = generated by two functions xi ^-nd X2- 
Taking xi and X2 to be homogeneous polynomials of degree 3 and 4 respectively, 
an elementary computations shows that 

H0T10T2 = Ho (3.1) 
+ {xuHo}+Hi (3.2) 

+ {X2,Ho} + ^{xiAxi.Ho}} + {xuHi} + H2 (3.3) 
+ h.o.t. (3.4) 



where the term (|3.2p is a homogeneous polynomial of degree 3, the term p.3p 
has degree 4 and h.o.t denotes higher order terms. We will construct a function 
Xi such that (|3.2p vanishes and we will show that there exists a X2 such that 
(|3.3p is reduced to normal form. 

To make precise the construction we need to split Hq as follows 

Ho = Hoo + Hoi, Hooiz):=Y,^^^ (3-5) 

fc 

„ Pl + ll 1 4sin^ {2WT2) ^ , 

^2 Wfc + 1 

Definition 3.1. A polynomial Z will be said to be in normal form if it Poisson 
commutes with ffoo, i-e. if 

{Hoo,Z} = 0. (3.6) 

Remark 3.1. In order to study the system with DEC we will always extend the 
system to a system defined for j ~ — (iV+ 1), ...,N with PBC, which is invariant 
under the involution qj i— > —q~j, Pj —p-j- The extension of the /3- model is 
obtained without modifying the equations, while the extension of the a- (3 model 
is given by the system 

9j Pj ' -Pj = -I3 - ai^il)] - (^Sjq^j - (3q^, (3.7) 
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where Sj is a discretization of the step function given by 



1 if j > 1 

,j .- ^^ if j = (3.8) 
-1 if.?<l 

and 

{^iq)j 2(7j - q^+i - (3.9) 

is the discrete Laplacian. The need of the introduction of the sequence Sj is at 
the origin of the finite smoothness of the solution in the DBC case. 



We are going to prove the foUowing 

Theorem 3.1. Assume a < 1/3. Then, there exists an analytic canonical 
transformation z = "^{ip), defined in a neighborhood of the origin 

HoT^ Ho{4') + Zii') + n{^) , (3.10) 

where Z is in normal form and the following holds true 

1) The remainder, the normal form and the canonical transformation are 
estimated by 

wxnms,. < c'Tzimi^, (3.11) 

WXzmis,. < CzMl, ' (3-12) 
\\^-rms,a < Criml, (3-13) 

i < s , cr > if PBC 



where 



s < I , (7 = if DBC 



and C-ji z,T ore constants independent of N . 
2) One has T = I + + TZt with xm given by 



and 



|2 



where the parameters vary in the range 



< si < s - i < 2, a = for DBC 
< si < s - i fsi < 2 , a>0 for PBC 



(3.14) 



!17^T(V')I1,_.,.,<C7^,M^MIV'IIL (3.16) 



(3.17) 
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3) The normal form has the following structure 

Z{tf},ij) = Zo + Zr 

where 

Zo:^lY.\^,\\ j:=l((3-faA (3.18) 

j 

and 

\\XzAn\s-s,,a<l^''Cr\ml. (3.19) 

with the parameters varying in the range (|3.17p . 



The proof of this theorem is divided into three parts. In the first one we 
will prove an abstract normal form theorem under the assumption that the non- 
linearity corresponding to the system (|2.4p is smooth. In the second part we 
will prove such a smoothness assumption. In the third part we will compute 
the first order term of the normal form and of the transformation and we will 
estimate the corresponding errors. 

We point out that there are 3 delicate points in the proof: the first one is 
to solve the homological equation (see lemma 13. 3|) ; the second one is to prove 
smoothness of the perturbation in the optimal space, and the third one is the 
actual computation of the main part of the normal form and of the canonical 
transformation. 



3.2 An abstract normal for theorem. 

First we recall that a homogeneous polynomial map F : Vs,a 'Ps,a of degree 
r is continuous and also analytic if and only if it is bounded, i.e. if there exists 
a constant C such that 

||F(z)||,^, <C||z||;^ yz&T,,,. (3.20) 

Definition 3.2. The best constant such that (|3.20[) holds is called the norm of 
and will be denoted by | F \s ,y One has 

I F I,,, := sup ||F(z)||,., . (3.21) 

Definition 3.3. A polynomial function /, homogeneous of degree r + 2, will be 
said to be of class Ti^ ^ if its Hamiltonian vector field ^/ is bounded as a map 
from Vs.cr to itself. 

Theorem 3.2. Let Hq be as above (see p.lip ) with a < 1/3. Assume that 
Hj S Ti-l ^ for some fixed s,a and for j ~ 1,2; then there exists an analytic 
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canonical transformation T = T10T2, defined in a neighborhood of the origin in 
Vs.a such that 

H oT = Ho + Z + n , (3.22) 

where Z is in normal form and the remainder, the normal form and the canon- 
ical transformation are bounded by 

\\Xniz)h,, < Cn\\z\\t,. , (3.23) 
11^2(^)11,,. < Cz\\z\\l, , (3.24) 
\\z-T{z)l^^ < Ct\\z\\1,, (3.25) 

with constants Cn^z,T depending only on a, | \s,cr, \ \s.a- 

The rest of this subsection wiU be occupied by the proof of theorem 13.21 
First we need some simple estimates. 

Lemma 3.1. Let f e TV,^^ and g e H2„, then {/,.g} G W+J^ and 

I X^f,g} |,,. <{r + r,+2)\Xf |,,. | X, (3.26) 

Proof. First remember that 

Xu,g} = [Xf.Xg] = dXfX, - dXgXj . (3.27) 

We recah now that, given a polynomial X of degree r + 1, there exists a unique 
(r + 1)— linear symmetric form X such that 

Xiz)=X{z,...,z) , 

then, (|3.27p is explicitly given by 

X{f .9} = i^ + ^)Xf{Xgiz), z, z, z) - in + l)Xg{Xf{z), z, z, z) (3.28) 

moreover from (|3.2ip one has 

X(zi,...,z.+i) < |^k.|kilL,...|k.+i|U (3-29) 

from which the thesis immediately follows. □ 

Remark 3.2. Let / G "HI „ then the corresponding vector field generates a flow 
in a neighborhood of the origin in 'Ps,cr- 

Lemma 3.2. Let x be of class Ti-l^, and let f G 'H'^^. Let $^ be the Lie 
transform generated by x, then each term of the Taylor expansion of f o $^ is 
a polynomial with bounded vector field. 
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Proof. Iterating the relation 

|(/°'i>*) = {x,/}°'i>* 

one gets that the Taylor expansion of / o cf)^ is given by 

l>0 

Then, the thesis follows from lemma [XT] □ 
A key role in the proof of theorem 13.21 is played by the so called homological 
equation, namely 

{Xo:Ho] + f,=Z, (3.30) 

where fj G TLl o- is a given polynomial, and Xj ^ T^i o-i ^ "^l a ^^'^ to be 
determined with the property that Zj is in normal form. 

Lemma 3.3. Consider the Homological equation p.30p with fj of class Til ^, 
with j = 1,2. Assume that a < then p.30p admits a solution Xj, Zj G Ti-i^^ 
with 

I ^x. < ^(Y^ I \s.- (3.31) 
Proof. First we rewrite the homological equation as 

(Lo + Li)x, =./,-Z, (3.32) 
where the operators Lq a-nd Li are defined by 

LaXj ■■= {Hm,X]} , LiXj := {-ffoi,Xj} 

and Hqo and Hqi are defined by (|3.5p . We will invert Lq and solve (|3.32p by 
Neumann series (see jBDGSOT] ). 

We begin by showing that the space HI j < 2 decompose into the sum of 
the kernel Ker(Lo) of Lq and of its range Im(Lo). Moreover, we show that Lq 
is invertible on its range. 

Given f GTii^ with j = 1,2 define 



Z ■■= ^£f{^\z))dt, 



1 '■^ 



(3.33) 

X := ^j^ [tf{^\z))-Z{^\z))]dt, (3.34) 

where 4'* is the flow of Xuoo a-iid T = 1 is its period. Then an explicit compu- 
tation shows that Z eKer(Lo), and that (see |BG93| ) 

Lox = / - ^ . (3.35) 
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Thus denoting by Q the projector on the kernel of Lq, and P = 1 — Q the 
projector on the range, one sees that the (|3.33p is a concrete definition of Q, 
while (|3.34p is the definition of L^^ restricted to Ini(Lo)- It remains to show 
that X € Ti-i Remark that, since is a canonical transformation one has 

Xz{z) ^ XQf{z)^^ j\^-*oXfO^S'){z)dt , (3.36) 

X^{z) EE X^-.p^.{z)^^ j\^-*oXpfO^S'){z)tdt , (3.37) 

From which it follows that 

I XQf |,., < I Xf |,,, , (3.38) 

where the last norm is the norm of L^^ as a linear operator acting on the space 
Hi,, and thus Qf G HI,, L^'Pf e HI,. 

We come now to the true homological equation p.32p . We look for a solution 
Xj = Pxj and Zj = QZj. Applying P or Q to p.32p . remarking that since 

[Lq, Li] = one has [P, Li] = [Q, Li] = 0, we get 

(io + Li)xj = P/, , Qfj = Z, . (3.40) 
The first of (|3.40p is formally solved by Neumann series, i.e. defining 

(Lo + Li)-i := ^(-1)" [L^^L^t L^' , X, ■= (^o + L.y'Pf, ■ (3.41) 

fc>0 

To show the convergence of the series in operator norm we need an estimate of 
To this end remark that, for any s,(t one has 

I Xhoi \sa < 2a , 

which using lemma [Ol implies j|ii|| < 2a{j + 1) < 6a. It follows that the series 
p.4ip converges provided a < 1/3, which is our assumption, and that 

||(^o + L,)-||<^^, (3.42) 

which concludes the proof. □ 
End of the proof of Theorem \3.S\ From lemma 13.31 one has that the solution xi 
of the homological equation with /i = Hi is well defined provided Hi G Hi , 
for some s,a. Then xi generates a Lie transform 71 which puts the system in 
normal form up to order 4. Then the part of degree four oi H o Ti takes the 
form 

/2 := \ {Xi, {xi.Ho}} + {xi.Hi} + H2^\ {xi.Hi} + H^ (3.43) 

which is of class Ti^ It follows that one can use the homological equation with 
such a known term and determine a X2 which generates the Lie transformation 
putting the system in normal form up to order 4. This concludes the proof of 
theorem 13.21 □ 
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3.3 Proof of the smoothness properties of the nonhnearity 

In this subsection we prove the following lemma 

Lemma 3.4. Let Hj, j = 1,2 be given by (|2.4[) . Consider the vector fields 
of the cubic and of the quartic terms of the Hamiltonian: they fulfill the 
estimates 

iiv / Ml / ^ II ii2 f ^ < cr > if PBC ,„ ... 

\\XHA^)h.<G^\\z\U^ , (i<.<|, a = ifDBC '^^'^^^ 

and 

\\XHA^)h^, < G2 \\z\\l^ , l<s , a>0 both cases (3.45) 

where we set Gj :~ | Xhj \s.cr- 

The proof will be split into two parts. First we show that it is possible to 
prove the result working on interpolating functions, and then we show that the 
"interpolating nonlinearities" have a smooth vector field when the parameters 
s, a vary in the considered range. 

Remark 3.3. Define T; as the map qj t-^ q'j^^ in the case of PBC, and [Ti((7)]j = 
q'^Sj and [T2{q)]j := qj in the case of DBC. Then the vector field Xif^ has only p 
components, moreover the norms are defined in terms of the Fourier variables, 
so we have to estimate the map constructed as follows 

— X! q\'"^^k{j)^Pk := ^/^pl 

3 

It is immediate to realize that the maps 1 and 2 are smooth (the frequencies 
are between 1 and 3) so it is enough to estimate the remaining maps. In turns 
the remaining maps essentially coincide with the map Ti(q) read in terms of 
standard Fourier variables (without the factors y^uJk)- These are the maps we 
will estimate. 

All along this section we will use a definition of the Fourier coefficients of a 
sequence not including the factors ^JuJk, namely we define qk by 

qj = ■ 

k 

We start by showing how to use the interpolation in order to make estimates. 
To this end we define an interpolation operator I by 

[/(<z)](x) :=^VMql.e^.(a:) (3.46) 

k 
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We also define a restriction operator R that to a function associates the 
corresponding sequence, by 

[Riu)], := ui^ij) (3.47) 

We remark that the operator R is defined on functions which do not necessarily 
have finitely many non- vanishing Fourier coefficients. 
Remark 3.4. With the definition (|2.14p one has 

ll^9ll,s,. = ll'z!l.,. . (3.48) 
Lemma 3.5. For any s > 1/2 there exists a constant Ce{s) such that one has 

\\Rul^„<Ce\\ul^^^ (3.49) 
Proof. Denote qj = {Ru)j] using the formula 



1 



efc(Mj) = efc+2(Ar+i)™(/^j) = 



(3.50) 



one gets 



N 



from which 

Let's define 
7rn,fc(s,cr) : = 



qk 



m 



N 



E 



Uk+2(N+l)m 



[k] 



k=-{N+l) 

(3.51) 

p'y\k\ 



[fc + 2(Ar + l)m]'^e'^l'=^+2(^+i)™l \[k + 2{N + l)m] J e'^l'=+2(^+i)'"l 



and replace p.5ip in the norm ll^ll^ g., then we get 



TV 

^ E 

fc=-(Af + l) 
TV 

k=-{N+l) 
N 



E [fc] 



2s2a\k\ 



k=^{N+l) 



E 
E 



Wfe+2(JV+1) 



7m,fc(g,0") , 
7m,/c(s, a) 



Uk+2{N+l)m 



< 



< E fE7L.(^,-))E 

fc=-(Ar+l) \m6Z / mSZ 



lc=-(W+l) 

< Ce{s)\\u\\ 



[^]2sg2^|fe| 

-2 7-— Tl"fc+2(Ar+l)m| 



< 



Indeed, since fc = —{N +1), . . . ,N, we can estimates the two factors of 'ym,k as 
follows: 
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|2(iV + l)m + fc| > |2(iV + l)|m| - \k\\ > 
which gives 



\k\, TO = 0, 

N + l>\kl TO^O, 



T|fc+2(Ar+l)m| 



< 1. 



W ^ I [2(JV+l)m] < J2^^ TO-0, 



which gives 



[k + 2{N+l)m]J - [2to]2«' 



< 



□ 



Corollary 3.3. Let T : M^iw+i) _^ r2(w+i) fte a polynomial map, assume that 
there exists an "interpolating polynomial map T'^ " such that T = RT'^I . If the 
map T'^ is bounded in some space H"''^ , with s > 1/2, then T is bounded in (.^ ^. 
Moreover one has 

I T I,,, <Cq\T' I,,, . (3.52) 

First wc define the interpolating maps we have to study. They are Tf{u) := 
in the case of PBC and Tl{u) = sgn(.T)u^(.T) and T!f(M) := in the case 
of DBC. Here we introduced the function 

( 1 if a; > 
sgn(a;) := < if a; = 
[-1 ifx<l 

The estimates p.44p in the case of PBC and p.45p for both boundary conditions 
are proved in Lemma lA.ll by a standard argument on the Sobolev norm of the 
product of two functions. We come to the estimate of Tf in the case of DBC. 

We win denote by the subspace of composed by the odd functions 
u{x) on [— TT, tt]. 

Lemma 3.6. For any 1/2 < s < 5/2, The operator 

T^{u) u'^sgn{x) 
is smooth from in itself and there exists C7{s) such that 

mui^<c\\\u\\i 
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Proof. We begin with the case 2 < s < |. First, observe that the function 
Tf (u) is odd when u is odd. we will prove the thesis by showing that the second 
weak derivative cP'T^{u) of T^{u) is in H^~^ . First remark that, by an explicit 
computation which exploits the fact that m(0) = one has 

d^T^{u) = 2[sgn(a;)(M'2) +sgn(a;)wd2u] . 

We show now that both terms are in H^^^ . The second term can be considered 
as the product of the function (Pu E H§^^ and of sgn(.T)u. This last function 
is of class , as it is seen by computing its derivative, namely 

d{sgn{x)u) = sgn(x)du{x) + S{x)u{x) — sgn{x)du{x) 

which clearly belongs to L^. From lemma lA.ll it follows that the product 

sgn.{x)u[d^u] e H''^'^. 

Concerning the term sgn(a;)(u'^), it can be considered as the product of 
^ ij^-i and of sgn(x), which is of class for all r < 1/2, as it can be 

seen by explicitly computing its Fourier coefhcients. Thus lemma [XT] gives the 

result. 

The case 1 < s < 2 is easier and works in a very similar way. Indeed, 
since u(x)sgn(a;) g and du{x) e iJ" with < s < 1, the derivative 
d{u'^ {x)sgn{x)) = 2u{x){du{x))sgii{x) belongs to iJ" with < s < 1, which 
gives the thesis. 

Quite different is the case ^ < s < 1, since by hypothesis no weak derivative 
exists for the function u{x). We exploit the following equivalent definition of 
the norm of the Sobolev space H^{[—tt, tt]) with real exponent s 



\u{x) - 

|^_y|l+2s 



■dxdy, 



(3.53) 



and the skew-symmetry of the periodic function u(x) G H^. We want to prove 
that g{x) := u(a;)sgn(a;) G H'^ with ^ < s < 1; more precisely 



\uix)sgn{x)\\l<4\\uix)\\l 



1 

-<.<1. 



(3.54) 



The symmetries of g{x) on the given domain allow to simplify the integral in 
(15351) 



= 2 



= 2 



< 



1.9(2;) -5(y)P 



Jo 



Jo 



Jo 



■dxdy 



\g{x)-g{y)\^ , \g{x) - g{y)\^^ 



y 



l+2s 



l+2s 
2 



\u{x) - u{y)\ 



\l+2s 



\u{x) - u{y)\'^ \u{x) - uiy)]"' 



\x + y 



l+2s 



\x-y 



l+2s 



■dxdy < 4 



dxdy = 
dxdy < 



□ 
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3.4 Computation of the normal form and of the transfor- 
mation. 

In this section we will concentrate on the case of DBC which is the most difficult 
one. 

Consider again the Hamiltonian (|2.4p , introduce complex variables £,j , r]j 
defined by 

"•-'^ <-^' 

and split Hq = Hqq + Hqi with 

-Hoo ^^^^EEi(^;r?),. , (3.56) 

i 



where (C;??)f2 := J2j ^jVj- 

Remark 3.5. The above splitting is different from the one introduced in p.Sp 
which had been used in the proof of lemma 13.31 and which was based on the 
Fourier variables. In particular one has {Hqq; Hqi} = 0, but {'Hoo;'Hoi} ^ 0. 

In the variables (^, 77) the flow of Hoo acts as follows 
The third order part of the Hamiltonian takes the form 



where Sj is the discrete step function defined in (|3.8p . The form of H2 will be 
given below. 
Denoting 

Co :-{Woo,.} , Ci :={Hoi;.} (3.60) 
we rewrite the homological equation for xi a-s follows 

iCo+Ci)xi=H, (3.61) 

which is solvable since the Kernel of £0 + Ci = Lo + Li on polynomials of third 
order is empty. The solution xi of p.6ip is unique and, as shown in Lemma 
13.31 exists. 

By a direct computation one has 

XI = C^'Hi - {Co + Cir^CiC^^Hi- (3.62) 

we are going to show that the second term is much smaller than the first one. 
Before starting, a couple of remarks are in order. 



3 
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Remark 3.6. The discrete Laplacian is f^-symmetric on periodic sequences 

< Aie,?7>^2=<e,Air/>^2 . (3.63) 

This is an immediate consequence of the fact that in Fourier variables it acts as 
a multipher by a real factor. 

Remark 3.7. In Fourier coordinates the discrete Laplacian Ai defined in p.9p 
acts as a multiplier by sin^ fc^. It follows that it has norm 1 when acting on 
anyone of the spaces Moreover, since 

|sin2(fc^)| < ri^"!, sie[0,2], fc^e[0,7r] 

one also has 

||AiCl|,_,,,. <M'^Mieil.,., 51 e [0,2]. (3.64) 
From ((3J4)) and ([S^T]) it follows 

\\Xno^{^,v)\\s-s,,o<C{a)fI''\mv)\Lo^ 0<,si<s-l. (3.65) 
Lemma 3.7. Assume a < ^, then xi = Xio + Xir with 

and there exists Cgia^Gi) such that 

Il^xi.(e,.)ll,,-.„.<C^8M^M1(C, (3.66) 

with 

< si < s - i < 2, cr = /or DBC 
^0 < si < s - i, si < 2 a>0 for PBC 

Proof. According to (j3.62p . let's define 
Xio :== Cq^Hi, xir ~{Cq + £i)"^£iXio = -(^o + ^O^^-CiXio- 

Since xw solves the homological equation {TioOiXio} = Hi^ it can be explicitly 
computed by 



Xio(0.) = f H,{^\tri))dt 
Jo 

a 
6V2 



which also implies 
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The thesis follows from ^^J^ and □ 

We move to the second Homological equation 

{X2,Ho} + H2 = Z (3.67) 

where = -^^2 + ^ {xii^i} can be split according to p.62p into = H20+H21 
with 

H20--^H2 + ^{Xw,H,}, H2i-^l{xir,H,}. (3.68) 
More explicitly, the leading term H2a is composed of 

^{Xio,i?i} = |jEfe' + 4ie^/, + 10C-??f-4iC,r,3+,,4)^ 



Before proceeding, it is useful to perform the change of variables 

= ^^^^ , tpj = E^fc^fc(-?) ■ 

which puts the quadratic part into diagonal form. This implies a modification 
of -^20 5 which however is of higher order and therefore will be included into the 
remainder terms. Indeed, the Lemma below shows that the difference between 
^ and ip is small 

Lemma 3.8. For any 2 > si > it holds true 

1 

S > -, cr > 0. 

Proof. By definition 

< 5:[fc]^^e^-nvi^.-i)^(l9\f + 1^^!^); 

k 

a Taylor expansion of the frequencies ujk gives 

' ' - 4 \2N + 2j - ^ 

which is the thesis. 

□ 
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Corollary 3.4. In terms of the variables tp^ip one has H2 = i?20 + H21 where 
2 

+ ^ E (^1 - + 6|^,f - 4i^, V^J + , (3.71) 

anii there exists 6*9(0, Gi) sitc/i t/iat 

\\XH,d^,^\^_^,^^^<C^^i''\ml„ (3.72) 

Just averaging (|3.7ip with respect to the flow $* it is now immediate to get 
the following Corollary. 

Corollary 3.5. The normal form Z is composed of two terms, Z — Zq + Zy, 
where the leading term Zq is smooth and reads 



3 

while the remainder is small 



II V / ; Ml ^ ^ SI II ; ii3 f < Si < s — i < 2 , (T = for DBC 

< C^9M- II^IL. |o <- ^ <s-^%,<2,a> Ofor PBC 

(3.74) 

Thus we have proved that the formula for Zq holds. The formula for xio 
implies that the canonical transformation has the structure p.lSp and this con- 
cludes the proof. 



4 Proof of Theorem 12.11 

To discuss this issue we first write the equations of motion of the first part of 
the normal form, namely of Hq + Zq, in the form 

i^A, = (AV'),-7V^,#jf , (4.1) 

where A is a linear operator which in the Fourier variables acts as a multiplier 
by Wfe = 1 + f /i2fc2 + Oin'^k'^), namely 

iM)k = ujki'k = (1 + ^/i'fc')V^fe + o(/) = (vWV')fc + o(m') ■ (4.2) 
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Take now an interpolating function for -0, in other words a function u such that 

ipj = eu{fij) , (4.3) 

where e is a small parameter representing the amplitude. Then, up to corrections 
of higher order, u fulfills the equation 

- iut ^ u + afi^Uxx - ie'^u\u\'^ , (4.4) 

which, up to a Gauge transformation and a scaling of the time introduced by 

u{x,t) = e^*ip{x,T) , T aj-P't, (4-5) 

gives the NLS equation 

\iPr = -<^xx + , 7:=-—. (4.6) 

a fjL'^ 

In order to get a hounded value of 7, from now on we take e = ^. 

We now compare an approximate solution constructed through NLS and the 
true solution of our Hamiltonian system. More explicitly, corresponding to a 
solution (^°(t) of the NLS with analytic initial datum, we define an approximate 
solution ■(/;" of the original model by 

iJ]{t) = ^le''^''{^lJ,a^lh) . (4.7) 

We also consider the true solution ip{t) of the Hamilton equation of the original 
model, with initial datum tp^ = fj.ip'^ {fij , 0) . 

From now on we will restrict to the case of DBC which is the complicate 
one. 

Wc first work in the variables in which H is reduced to the normal form 

H = Ho + Zo + Zr+TZ. (4.8) 

Lemma 4.1. Let tjj be the solution of the equations of motion of (|4.8p with 
initial datum ip{0) = T~^(ip^) and let ip"" be as defined in (j4.7p . then 

with 

Proof. Observe that the NLS equation for tp'^ may be rewritten as 

r = xh.,sA^-) + ^zir) + XTzir) - xnAr) (4.9) 

where 

7^l ~ {Ha-HNLSfl) + Zr+n 



.2,2 



2 
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so that Xti^ fulfills the esthiiate 



< 



\\xndr)\\.,o < c^i^'\\r\\s+4,o + Cr^^''\\r\\:+suo + ^'T^\n.o 

We compare -0" with the full solution ip of the equation 

with initial datum ?/'o = T~^(?/''^), whose difference from i/j'^ (initial datum for 
tp"') is controlled by 

So we apply the Gronwall lemma (see lemma |A.3[) with 

A := Xho , P = Xz + X-ji , R ~ X-jz^ , 
obtaining that the error ipi := ip — ip*^ from the NLS dynamics satisfies 

i;i=AtPi + [P{r + i'i)-P{in]+XnAi^") (4.10) 
and is estimated by 

||V'i(i)L.,<Cr/.'e^«^^* + ^(e^«^^*-l) <C7A.^ \t\<T/^,^ (4.11) 

where Cq := 6Cz- 

□ 

When we go back to the original variables, the solution z may be split as 
z = r (V-) = T(V''' + V'l ) = M^" + M^zi + /i^zz 
where we have defined 

a r ^l+^X0l(V'°) Z - //Z" - 

z := — , zi := ^ , Z2 :— ^ . i^.i^j 

More precisely, we claim that it holds 
Lemma 4.2. We have 

lkllL,0<Cl, ||^2|L,o<C^2, 

MP to times s.t. \t\ < ^ . 
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Proof. The first inequality comes directly from the lemma H?T1 Concerning the 
second one, we remark that 

z = r(V'" + ^i) = ^" + v^i + [r(r + ^i)-(V'" + ^i)] = 

= r + v^i + + v^i) + o{\\r + Mlo) = 

= r + V'l + x^^ir) + o(ii^'iii.,o \\r\Lo + \\r + mIo) 

(by differentiability of X^-^ and Lagrange mean value theorem). Finally from 
(|3.66p with si = 1 we have 

□ 

We now analyze the first correction zi. To this end we analyze separately 
its two terms. 

First remark that, from (|4.10[1 one has ijji = ipio +0{\t\ H^'^Hg o)' where V'lo 
solves the equation ipio = AipiQ with initial datum T^^{ip'^) — 'ip^. Thus we have 



(4.13) 



We now analyze the other term. To this end, with the aim of considering the 
short time dynamics, we rewrite the equation (|4.9|1 as 



^'^ = zV" + 0(m') ^ V;'^ = e'V° + 0(M'|t|) 
Thus exploiting the differentiability of ^^oi have 

M^zi = -c^*X^o, (^°) + X^o, (c'V°) + Oif,'\t\) . 

Since 



(4.14) 



e^*X^„(V'°) = - 



6V2 

a 



3e^V'o+6ie^'lV'o|' + e^^^o 



At- 



3V'^e^'* + 6i|Vo| 



-2ii 



which yields to 
a 



6V2 



3(e 



At „2itN„;,2 



)^o' + 6i(e^ -l)l^o|' + (e 



In the case of tJjq — i/izo (zero velocity initial datum) we have 



6\/2 



f4e^* - 3e2i* _ g-si* _ 6ie 



At 



6i] z2 + 0(/i|t|) 



which gives immediately the thesis. 
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A Appendix: a few technical lemmas. 



Lemma A.l. Let u e H'^'^ and v S i?*-'^ with s > i and s>r>0, a>0. 
Then there exists C — C(r, s) such that the following inequality holds 

\\uv\l^„<C\\u\l,J\vl^^ (A.l) 

Proof. In this proof it is useful to use the expansion of u and v on the complex 
exponentially. Thus we will write u{x) = X^fcGZ ^'^'^ remark that 



if in the definition of the norm cf (|2.14p we substitute such coefficients to the 
coefficients on the real Fourier basis, nothing changes. This is due to the fact 
that both the basis of the complex exponentials and the real Fourier basis are 
orthonormal. The advantage is that in terms of the complex exponentials the 
product is mapped into the convolution of the Fourier coefficients, thus we have 
simply to estimate the norm of the function whose Fourier coefficients are 

{u *v)k U]-kVk ■ (A. 2) 

'j 

As a preliminary fact we define the quantities 

and prove that there exists a constant C(s, r) such that 

E i < ^ (A.3) 

To obtain (|X3|) we need some preliminary inequalities. For any positive a and 
b one has 

(a + 6)2'- < 22'- max {a; b}'^'' < r{a'^'' + b^'') (A.4) 
and for any j and fc in Z 

1 1 1 1 .A r^ 



[j -fc]2'-[fc]2'i-2r (min{[j-fc],[fc]})2'' [j-kf' [kf 
From (|A.4p . (jA.Sp and [j] < [j - fc] + [fc] it follows 



< 4'' ( ^— + 1 < 

W [j-fc]2'-[fc]2«-2r 



< 4"^ 

which gives (|A.3p with 



1 1 1 



[fc]2'' [j-fc]2« [fc]2« 
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Hence 



i k 



< 



which concludes the proof. □ 



Wc state here a version of the GronwaU Lemma which is suited for our 
estimates. First we recall the following lemma. 

Lemma A. 2. Let x : [0,T] V be a differ entiahle function and V a Banach 
space. Assume that Wt g [0, T] it fulfills the integral inequality 

\\x{t)\\<K+[ (a||a:(s)|| +6)ds (A.6) 
Jo 

with a, b real and non negative parameters, then 

\\x{t)\\ < e-^'A' + - (e"* - 1) . (A.7) 
a 

The lemma we use in sect. |4]is the following one. 
Lemma A. 3. Let z{t),z°-{t) gV, t G [—T,T] be respectively the solutions of 

jz = Az + P{z), j z" ^ Az'' + Piz") - Riz"), 

\z{0) = zo ' \z-{0) = zS 

where A is the generator of a unitary group in V , and \\z°'{t)\\ < C. Assume 
also that the non-linearity P has a zero of third order at the origin and that for 
all t € [-T,T] and all z with \\z\\ < 2C 

\\Piz)\\<p,\\z\\\ \\dPiz)\\<3p,\\z\\\ 

and the remainder R is estimated by 

WRiz'^mi < P2, yte[-T,T]. 

Let 5 z — z"" , then the following estimate holds 

II^WII < ||<5(0)||e''^* + g(e''^*-l), tG[-T,T]. (A.8) 
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Proof. The difference S{t) is solution of the differential equation 
S^AS+ +S)- + Riz^it)) , 

by Duhamel formula one has 

S{t) = e^*<5(0) + e^* /* e"^^ [(Piz^is) +5)- P(z"(s))) + i?(^"(s))] ds 

^0 

Using Lagrange mean value theorem to estimate P{z°'{s) + S) — P(z°(s)) and 
the fact that A is unitary one has 

II^WII < 11-5(0)11+ f\3p,\\S{s)\\+p2]ds, 
Jo 

which fulfills (|A.6|) . from which the thesis follows. □ 
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